Reply to the comment of Kolakowska and Novotny 



o 
o 

O 
CD 

Q 



Claudio Horowitz and Ezequiel V. Albano 
Institute de Investigaciones Fisicoqufmicas Teoricas y Aplicadas (INIFTA), Universidad Nacional de 
La Plata, CONICET, CCT-La Plata CONICET; Sucursal 4, CC 16 (1900) La Plata, Argentina 

(Dated: July 21, 2011) 

In an early paper (Horowitz and Albano, Phys. Rev. E., 73 031111 (2006)) we studied growing 
models, generically called X/RD, such that a particle is attached to the aggregate with probability 
p following the mechanisms of a generic model X, and at random [Random Deposition (RD)] with 
probability (1 — p). We also formulated scaling relationships that are expected to hold in the limits 
p — > and L — > oo, where L is the sample side. In the previous comment, Kolakowska and Novotny 
(KN) stated that our scaling hypotheses does not hold. Here, we show that the criticisms of KN are 
outlined by analyzing data out of the proper scaling regime and consequently they are groundless 
and can be disregarded. 

PACS numbers: 89.75.Da, 81.15.Aa, 68.35.Ct 
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In the previous paper Kolakowska and Novotny (KN) 
[l[ comment our early paper on the " Dynamic proper- 
ties in a family of competitive growing models" [2[. Our 
paper, as well as a series of previous manuscripts on a 
related issue @, 0, @| , addresses the properties of a wide 
variety of growing models, generically called X/RD, in- 
volving the deposition of particles according to compet- 
itive processes, such that a particle is attached to the 
aggregate with probability p following the mechanisms 
of a generic model X that provides the correlations, and 
at random [Random Deposition (RD)] with probability 
(1 — p). A related study on that topic has also been 
published by Braunstein et al. [6J. 

The comments of KN are based essentially on two 
statements, both of them erroneously attributed to us, 
which read as follows: 

Sla) "The claim is made that at saturation the surface 
width w(p) obeys a power-law scaling w(p) oc l/p s , where 
5 is only either S = 1 or 5 = 1/2, which is illustrated by 
the models where X is ballistic deposition and where X 
is RD with surface relaxation." (Taken from the Abstract 
of the comment of KN [l| ) . A slightly different version of 
this statement can also be found in the comment of KN 
[l|, few lines below equation (1), namely: 

Sib) "The new claim that is being made in Ref.[l] (i.e. 
reference 0) is that a nonuniversal and model- dependent 
exponent S in Eq.(l) must be only of two values, either 
S = 1 or S = 1/2, for models studied in Ref. [1] (i.e. 
reference [2j)." 

On the other hand, the second statement, taken from 
the Abstract of the comment of KN [l[, reads: 

S2) "Another claim is that in the limit p — > 0, for 
any lattice size L, the time evolution of w(t) generally 
obeys the scaling w(p) ex (L a /p s )F(p 2S t/L z ), where F is 
Family- Vicsek universal scaling function." 

Concerning both Sla) and Sib), our answer is that 
the statement is taken out of context from our paper, so 
it is incomplete and leads the reader to confusion. In 
fact, in our paper we explicitly state in many places that 



the values 6 = 1 or 6 = 1/2 only "hold" in the p — ► 
limit. In fact, in our paper we state in the Abstract, 
the paragraph before the title of Section V, in the title of 
Section V, and in the Conclusions that the universality in 
the 6 and y exponents should "hold" in the p — > limit. 
Also, we have theoretically found those values of 6 in the 
"p — > 0" limit by using a correspondence between two 
neighboring sites in the discrete model ((h(i) — h(i + 1))) 
and two types of random walks. 

Within this context it is worth mentioning that Braun- 
stein et al. @ have determined exactly that S = 1/2 for 
a competitive model between Ballistic Deposition (BD) 
and RD, while they found 6 = 1 for the competition 
between Random Deposition with Surface Relaxation 
(RDSR) and RD. Of course, Braunstein et al. [6j clearly 
state that the "scaling exponents derived are exact for 
p — » 0." They also recall that "at finite p, we expect devi- 
ation from the exact scalings which indeed were observed 
numerically (in reference [4|." See the last paragraph of 
Section II, page 2. 

Figures 1 (a) and 1 (b) of the comment of KN [lj show 
the behavior of various models, with X = RDSR, X = 
BD and some variants. Results are shown for a wide 
range of p, actually most results correspond to p > 0.3, 
i.e., far away from the correct scaling regime given by the 
p — > limit. Subsequently, by describing the figures, 
they state that "in special cases an approximate power- 
law w{p) oc 1/p may be observed, however, this is not a 
principle." From our point of view, it is obvious that one 
would not expect a nice power-law fit of the data within 
the whole range of p, but "only in the limit of p — > 0." 
Also, a careful inspection of figures 1(a) and 1(b) of the 
comment nicely shows that the power-law predicted in 
our paper fits very well the data for p — » 0!. So, the data 
shown in figure 1 of the comment of KN [l[ seem to be 
correct, but they are far away from the right scaling limit. 
The figure clearly suggests that for data taken correctly, 
namely, for p << 0.1 as stated in our paper [2J, they 
would certainly obey the proposed scaling behavior. 
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In a related context, also at the end of the subsection 
Saturation, KN state that [l[ "The other two examples 
shown in fig. 1 defy a linear fit. In these cases there is no 
power law of the type claimed in Ref.[l]" (i.e., reference 
[2j ) . This absence of power-law scaling in p is also evident 
in fig. 4 of Ref.[l] (i.e. reference Q). We notice that in 
contrast to that opinion, an excellent power-law behavior 
can be observed in figure 4 of our paper [2| , but of course 
in the right scaling regime, namely, for p — > 0. 

Concerning S2, our answer is also that the statement 
is incomplete and incorrectly formulated, so it leads the 
reader to confusion. In fact, in reference Q we clearly 
state that S2 only holds, as is well known for the stan- 
dard scaling relations, in the L — * oo limit, but not "for 
any lattice size L" as in the misleading statement of KN. 
The right scaling regime for the sample size is stated in 
many places of our paper, e.g. along the explanations of 
Section II, as well as the explicit scaling indications writ- 
ten in equations (6), (7) and (8) Q. We only mentioned 
in Section III that there is numerical evidence of negli- 
gible finite-size corrections to the value of the exponent 
5. Obviously, this property of 5 is not equivalent to the 
statement formulated by KN in a general way. 

All discussed issues and the remaining topics of the 
comment suggest that KN did not understand the con- 
cept of scaling in the p — > and L — > oo limits, and 
therefore they are unable to understand either the correct 
content of our paper or their own results. For example, 
in the subsection The RD limit, another comment of 
KN reads: "Another claim of Ref.[l] is that Eq.(l) with 
the power-law prefactors p s where (5 = 1 or 1/2) would 
prevail in the limit p — > 0, and that such a scaling would 
be universal. We tested these claims in simulations of 
RD+BD models and found the evidence to the contrary 
(fig. 2-3)." Again, figures 2 and 3 of the comment of KN 



show the behavior of RD+BD models using p in the range 
[0.1, 1] !, which of course is not at all close to the limit 
p -> 0. 

We notice that it is easy to understand why KN claim 
that our scaling did not fit the data very well. This is 
the obvious result that one may always obtain just by 
working with data taken out of the scaling regime but 
using values for the exponents corresponding to the cor- 
rect regime. In fact, as we have already shown in our pa- 
per, as well as in previous work [3j, jjj, when using small 
samples, data collapse can only be obtained by using the 
effective values of the exponents according to the range 
of p and L used. 

In summary, our paper shows that for some models 
and within the p — > and L — > oo limits, the proposed 
scaling is a universal principle. On the other hand, in 
their comment, KN have shown that data taken out of 
the correct scaling limit may depart from scaling. 

We acknowledge financial support from the Argen- 
tinian Science Agencies CONICET and ANPCyT, and 
the UNLP. 
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